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We show that the coherence of an electron spin interacting with a bath of nuclear spins can exhibit 
a well-defined purely exponential decay for special ('narrowed') bath initial conditions in the presence 
of a strong applied magnetic field. This is in contrast to the typical case, where spin-bath dynamics 
have been investigated in the non-Markovian limit, giving super-exponential or power-law decay of 
correlation functions. We calculate the relevant decoherence time T2 explicitly for free-induction 
decay and find a simple expression with dependence on bath polarization, magnetic field, the shape 
of the electron wave function, dimensionality, total nuclear spin /, and isotopic concentration for 
experimentally relevant heteronuclear spin systems. 
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I. INTRODUCTION 

There are many proposals to use the spin states of elec- 
trons in confined structures for coherent manipulation, 
leading to applications in quantum information process- 
ing and ultimately, quantum computation.— A se- 
ries of recent experiments on such spin states in quantum 
dots,— 1£ electrons bound to phosphorus donors in silicon, 8 
NV centers in diamond ) 4 ' 9 ' 10 and molecular magnet o 11 ' 12 
have shown that the hyperfine interaction between con- 
fined electron spins and nuclear spins in the surrounding 
material is the major obstacle to maintaining coherence 
in these systems. 

Previous studies of this decoherence mechanism have 
pointed to the non-Markovian nature of a slow nuclear- 
spin environment, leading to non-exponential coherence 
decay, 1 ?' 14 ' 1 ^ 16 ^ 7 ' 18 ' 19 ' 20 ^ 1 - 22 ' 23 - 24 ^ 5 These results sug- 
gest that it may be necessary to revise quantum error cor- 
rection protocols to accommodate such a 'nonstandard', 
but ubiquitous environment.— In the present work, we 
show that virtual flip-flops between electron and nuclear 
spins can lead to a well-defined Markovian dynamics, giv- 
ing simple exponential decay in a large Zeeman field and 
for particular initial conditions (a 'narrowed— nuclear- 
spin state). Moreover, we calculate the decoherence time 
T2, revealing the dependence on many external parame- 
ters for a general system. 

The rest of this paper is organized as follows: In Sec. 
HIl we introduce the Hamiltonian for the Fermi contact 
hyperfine interaction and derive an effective Hamiltonian 
for electron spin dynamics which is valid in a strong mag- 
netic field. In Sec. IHII we present the Markov approxima- 
tion and its range of validity, giving an analytical expres- 
sion for the decoherence time T2 . We also give bounds for 
the non-Markovian corrections to our expression. Sec. IIVI 
gives a discussion of the decoherence rate for a homonu- 
clear system and in Sec. |V]we generalize these results for 
a heteronuclear spin bath, providing explicit analytical 
expressions for T2 within our Born-Markov approxima- 
tion. We conclude in Sec. IVII and present additional 



technical details in Appendices [ATP] 

II. HAMILTONIAN 

We begin from the Hamiltonian for the Fermi contact 
hyperfine interaction between a localized spin- 1/2 S and 
an environment of nuclear spins, 

H M = bS z + 6^ 7fc 4 2 + S-h; h = ^A fc Ife. (1) 

k k 

Here, Ik is the nuclear spin operator for the spin at 
site k with associated hyperfine coupling constant Ak, 
b = g* /j,bB is the electron Zeeman splitting in an applied 
magnetic field B and 7^ is the nuclear gyromagnetic ra- 
tio in units of the electron gyromagnetic ratio (we set 
h = 1): 7fc = gi k ^N I ' 9* Ms- F° r an electron with en- 
velope wave function ^(r), we have Ak = voA lh \ip(rk)\ 2 , 
where A lk is the total coupling constant to a nuclear spin 
of species ik at site k and vq is the volume of a unit 
cell containing one nucleus. For convenience, we define 

Vi is the relative concentration 
of isotope i. The envelope function ip(r) of the bound 
electron has finite extent, and consequently there will be 
a finite number ~ N of nuclei with appreciable A^. For 
typical quantum dots, N ~ 10 4 — 10 6 , and for donor im- 
purities or molecular magnets, N ~ 10 2 — 10 3 . In Eq. 
(JTJ) we have neglected the anisotropic hyperfine interac- 
tion, dipole-dipole interaction between nuclear spins, and 
nuclear quadrupolar splitting, which may be present for 
nuclear spin I > 1/2. The anisotropic hyperfine interac- 
tion gives a small correction for electrons in a primarily 
s-type conduction band^& such as in III-V semiconduc- 
tors or Si. Nuclear dipole-dipole coupling can give rise 
to dynamics in the spin bath, which can lead to electron- 
spin decay due to spectral diffusion on a time scale found 
to be T M - 10 - 100 (is for GaAs quantum dotsi* 4 - 2 ^ 
These times are one to two orders of magnitude longer 
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than the T 2 we predict for a GaAs quantum dot carrying 
N = 10 5 nuclei (see Fig. [31 below). For smaller systems, 
we expect the decay mechanism discussed here to dom- 
inate dipole-dipole effects substantially. The quadrupo- 
lar splitting has also been measured for nanostructures 
in GaAs, giving inverse coupling strengths on the or- 
der of 100 /zs^i comparable to the dipole-dipole coupling 
strength, so quadrupolar effects should become relevant 
on comparable time scales. 

For large b, we divide Hm — Hq + Vg into an un- 
perturbed part H that preserves S z and a term Vg — 
\ (S+h- + S-h + ) that leads to energy non-conserving 
flip-flops between electron and nuclear spins.— We elim- 
inate Vg to leading order by performing a Schrieffer- 
Wolff-like transformation: H = e s Hhfe~ s « H = 
Hq + t; [S, Vg], where S — -^Vg, and l_o is the unperturbed 
Liouvillian, defined by LoO — [Hq,0]. The resulting ef- 

(see Appendix |B|) 



H = (lo + X) S z 



D. 



(2) 



The operators to and D are diagonal with respect 
to a product-state basis of I|-eigenstates (^) fc \Ii k m,k), 
whereas the term X is purely off-diagonal in this ba- 
sis, leading to correlations between different nuclei. We 
neglect corrections to the diagonal part of H of order 
~ A 2 /Nb, but retain the term of this size in the off- 
diagonal part X. This approximation is justified since, 
as we will show, the bath correlation time t c is much 
shorter than the time scale where these diagonal correc- 
tions become relevant for sufficiently large Zeeman split- 
ting b 3> A, where a Born-Markov approximation is valid: 
t c ~ N/A -C Nb/A 2 . In addition, we ignore corrections 
to X that are smaller by the factors Ak/b ~ A/Nb <C 1 
and 7^ ~ 10 -3 . Under these approximations, the various 
terms in Eq. ([2|) are given by (see also Appendix |B} : 



X 



b + h z , D~bJ2lkI z k , 

k 

1 ^ A ^ A l 



2 - u 



(3) 
(4) 



of the transverse electron spin (S+) t are described by the 
exact equation of motion^ 5 - 

(S+) t = iu n (S+)t - 1 I df E(t - (S+) t , , (5) 
Jo 

(6) 



E(t) = -iTiS + Le- lWLt QL\n) (n\S- 

Here, L and Q are superoperators, defined by their ac- 
tion on an arbitrary operator O: LO — [H,0], QO = 
(1 — |n) (n\ Ttj) O, where Trj indicates a partial trace 
over the nuclear spin system. 

To remove fast oscillations in (S+) t we transform 
to a rotating frame, in which we define the coher- 
ence factor xt = 2exp [— i(u) n + Aw)t] (S+) t and asso- 
ciated memory kernel £(t) = exp [— i{u n + Aw)t] £(t), 
with frequency shift determined sclf-consistently through 
Aw = —Re f °° dtT,(t). Additionally, we change integra- 
tion variables to t = t — t' . The equation of motion for 
xt then reads 



x t 



i 1 drl](r)xt_ 
Jo 



(7) 



If £(t) decays to zero sufficiently quickly^ 4 - on the time 
scale t c <^T 2 , where T 2 is the decay time of x t , we can 
approximate x t - T ~ %t and extend the upper limit on 
the integral to t — > 00 (Markov approximation), giving 
an exponential coherence decay with a small error e(t): 



x t = exp {-t/T 2 ) x Q + e(t), 



1 

T 2 



-Im / dfE(t). (8) 
'0 



The non-Markovian correction e(t) can be bounded 
precisely if £(t) is known) 3 ^ 



dt' 



dt"t(t") 



(9) 



Eq. @ gives a hard bound on the validity of the Markov 
approximation, and consequently, any corrections to the 
exponential decay formula. Fig. [1] demonstrates an ap- 
plication of Eqs. ijHJ) and ^ for decay in a homonuclear 
spin system, which we discuss below. 



IV. HOMONUCLEAR SYSTEM 



III. MARKOV APPROXIMATION 

For large 6, H^g leads only to incomplete decay of the 
longitudinal spin {S z ) t ^ However, it is still possible for 
the transverse spin {S+) t to decay fullyJ-^ through a pure 
dephasing process, which we now describe in detail. We 
assume that the electron and nuclear systems are initially 
unentangled with each other and that the nuclear spin 
system is prepared in a narrowed state (an eigenstate of 
the operator u>: to\n) = uj n \n)) through a sequence of 
weak measurements ; 27 ! 33 ' 34 polarization pumping, 35 fre- 
quency focusing under pulsed optical excitation,- 3 ^ or by 
any other means. For these initial conditions, dynamics 



If only one spin-carrying nuclear isotope is present, 
7fc = 7, independent of the nuclear site. We then 
approximate £(t) to leading order in the perturbation 
V — XS Z (Born approximation, see Appendix [C]) by ex- 
panding Eq. ([5]) through iteration of the Dyson identity: 
e - i LQt = e - l L Qt _ i f^t'e-f-oQCi-tOLyQe-^Q*', where 
LyO = [V,0]. Higher-order corrections to the Born ap- 
proximation will be suppressed by the small parameter 
A/u! n J^ Inserting the result into Eq. ((5J) we find: 

— = Re / dte- tAu]t (X(t)X) ; X(t) = e - wt Xe tuJt . 
T2 Jo 

(10) 
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Figure 1: (Color online) Exponential decay xt = exp( — t/T-z) 
(solid line) and maximum error bounds xt ± |e(i)| max (dashed 
lines), found by numerical integration of Eq. © with param- 
eters for a two-dimensional quantum dot (before Eq. (|14|l '). 
/ = 3/2 and A/b = 1/20. For comparison, we show the de- 
cay curves for super-exponential forms exp [— (t/T2) 2 ] and 
exp [— (t/T2) 4 ] (dot-dashed lines) and rapidly decaying bath 
correlation function C(t)/C(0) (dotted line, see Eqs. (|10[1 and 

(TO). 



Here, (•■■) = (n\ ■ ■ ■ \n) denotes an expectation value 
with respect to the initial nuclear state. Eq. (TIT)]) re- 
sembles the standard result for pure dephasing in a weak 
coupling expansion, where X(t) would represent the bath 
operator in the interaction picture with an independent 
bath Hamiltonian. However, for the spin bath there is 
no such weak coupling expansion, and X(t) appears in 
the interaction picture with u>, the same operator that 
provides an effective level splitting for the system. Ad- 
ditionally, the general result for a heteronuclear system 
including inter-species flip-flops cannot be written in such 
a compact form^ 

Previously, it has been shown that a Born-Markov ap- 
proximation to second order in Vff leads to no decay*^ 
In contrast, a Born-Markov approximation applied to the 
effective Hamiltonian leads directly to a result that is 
fourth order in Vs [Eq. (fTOf ]. describing dynamics that 
become important at times longer than the second-order 
result. It is not a priori obvious that the effective Hamil- 
tonian, evaluated only to second order in Vg, can be used 
to accurately calculate rates to fourth order in Vs- We 
have, however, verified that all results we present here 
are equivalent to a direct calculation expanded to fourth 
order in Vg at leading order in A/b <C 1<2& 

If the initial nuclear polarization is smooth on the scale 
of the electron wave function, the matrix elements of op- 
erators like I k 1^ can be replaced by average values. Ne- 
glecting corrections that are small in A/Nb <C 1, this 
gives (see also Appendix |D|) : 



C(t) = (X(t)X(0)) 



-A,)t 



(11) 



Above, we have introduced the coefficients c± = 1(1 + 
1) — ((m(m ± 1))) and the double angle bracket indicates 
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Figure 2: (Color online) Geometrical factor f(d/q) from Eq. 
(|13p . where d — 1,2,3 is the dimension and q characterizes 
the electron envelope function ?/<(r) = tp(0) exp [— (r/ro) 9 /2]. 

an average over 1^ eigenvalues 

In the limit N ^ 1 we can include the term k = I in 
Eq. and perform the continuum limit Efe — > J dk 

with small corrections. For an isotropic electron wave 
function of the form ip(r) = ip(0)e^i r ^ r °) ^ 2 contain- 
ing N nuclei within radius ro in d dimensions, the hy- 
perfine coupling constants are distributed according to 



A k = A exp 



(k/N) 



q/d 



where k is a non-negative 
according to 



index, and we choose Aq to normalize A^ 
A = / °° dkAifS- (see also Appendix [AJ ■ 

After performing the continuum limit, C(t) will decay, 
with characteristic time r c given by the inverse band- 
width of nuclear flip-flop excitations r c ~ 1/Aq ~ N/A. 
For large b, I/T2 will be suppressed due to the small- 
ness of X (see Eq. Q), whereas r c remains fixed. At 
sufficiently large b, it will therefore be possible to reach 
the Markovian regime, where r c is short compared to T2: 
t c /T 2 1. Evaluating the time integral in Eq. (fTU)) . 
we find the general result to leading order in A/uj n (see 
Appendix [D]) : 



1 

fir) - 



7T 



c+c-f 



4 

r { 3 



i\ 2r - 1 r(2 



A_y a 

- - 1 



[r(r) 



r > 1/2. 



(12) 
(13) 



In Eq. (Tj"2")) . A/N sets the scale for the maximum decay 
rate in the perturbative regime, the coefficients c± set 
the dependence on the initial nuclear polarization p (e.g., 
with / = 1/2, we have c + c_ = (1 -p 2 )/4), A/u n < 1 
gives the small parameter which controls the Born ap- 
proximation, and f(d/q) is a geometrical factor (plotted 
in Fig. [2]). f(d/q) is exponentially suppressed for d/q > 1 
(f(r) (x (l/S) 2r - l (l/r) r , r > 1), but f(d/q) -> 00 for 
d/q — 1/2 — > + . Due to this divergence, no Markov 
approximation is possible (within the Born approxima- 
tion) for d/q < 1/2. We understand the divergence in 
f(d/q) explicitly from the asymptotic dependence of C(t) 
at long times: C(t) cx l/t 2d /i, t > N/A, d/q < 2. 44 
Surprisingly, there is a difference of nearly two orders of 
magnitude in I/T2 going from a two-dimensional (2D) 
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find (see also Appendix [Dj) : 
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Figure 3: (Color online) Decay rates for an In^Gai-^As 
quantum dot with In doping x = 0.05. Here, we have as- 
sumed N = I0 5 and used values of Vi and A 1 for GaAs 
from Ref. gj: A ?5As = 86 fieV, J 4 6£ ' Ga = 74/ieV, A ?lGa = 
96 fieV, u-5 As = 0.5, ^69 Ga = 0.3(1 — x), U7i Ga = 0.2(1 — a;). 
The hyperfine coupling for In in InAs was taken from Ref. 
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quantum dot with Gaussian envelope function (d = 2, 
q = 2, d/q=l) to a donor impurity with a hydrogen-like 
exponential wave function (d = 3, q = 1, d/q = 3), if all 
other parameters are fixed (see Fig. [2]). 

We now specialize to an initial uniform unpolarized 
spin bath, which is nevertheless narrowed: oj\n) = b \n), 
with equal populations of all nuclear Zccman levels (i.e., 
(( m )) = and ((m 2 )) = + 1)). For a 2D quantum 



dot with a Gaussian envelope function (d = q 
find, from Eqs. (12|) and US]): 



i _ 7T n(i + i)Ay a 

T 2 ~ 3 I 3b ) N' 



2) we 



(14) 



There are two remarkable features of this surprisingly 
simple result. First, the condition for the validity of 
the Markov approximation, T 2 > t c ~ N/A will be 
satisfied whenever A/b < 1, which is the same con- 
dition that validates a Born approximation. Second, 
1/T 2 has a very strong dependence on the nuclear spin 
(I/T2 oc I 4 ). Thus, systems with large-spin nuclei such 
as In (/in = 9/2) will show relatively significantly faster 
decay (see, e.g., Fig. [3]). 



T< = — = vf — 

XI? 1 3 



X(X + 1)A' 
36 



A* 



TP (15) 



The quadratic dependence on isotopic concentration 
is particularly striking. Due to this dependence, electron 
spins in GaAs, where Ga has two naturally occurring 
isotopic species, whereas As has only one, will show a 
decay predominantly due to flip-flops between As spins. 
This is in spite of the fact that all isotopes in GaAs have 
the same nuclear spin and nominally similar hyperfine 
coupling constants (see Fig. [3]). Interestingly, we note 
that the relatively large flip-flop rates for In and As, due 
to large nuclear spin and isotopic concentration, respec- 
tively, may partly explain why only Ga (and not In or 
As) spins have been seen to contribute to coherent effects 
in experiments on electron transport through (In/Ga)As 
quantum dots.— The same effect may also explain why 
polarization appears to be transferred more efficiently 
from electrons to As (rather than Ga) in GaAs quantum 
dots42 



VI. CONCLUSIONS 

We have shown that a single electron spin can ex- 
hibit a purely exponential decay for narrowed nuclear- 
spin bath initial conditions and in the presence of a suffi- 
ciently large electron Zeeman splitting b. This work may 
be important for implementing existing quantum error 
correction schemes, which typically assume exponential 
decay of correlation functions due to a Markovian envi- 
ronment. In the limit of large Zeeman splitting b > A, 
where a Born-Markov approximation is valid, we have 
found explicit analytical expressions for the decoherence 
time T2, giving explicit dependences on the electron wave 
function, magnetic field, bath polarization, nuclear spin, 
and isotopic abundance for a general nuclear spin bath. 
Moreover, within the Born-Markov approximation, we 
have found a divergence in the decoherence rate 1 /T 2 for 
a one-dimensional quantum dot, indicating a breakdown 
of the Markov approximation in this case. 
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V. HETERONUCLEAR SYSTEM 

For sufficiently large b (|'yfe — Tfe' I ^ > \A k - A k >\ ~ 
A/N), heteronuclear flip-flops between two isotopic 
species with different 7fc are forbidden due to energy con- 
servation. In this case, 1 /T 2 is given in terms of an inco- 
herent sum, I/T2 = r = J2i where Tj is the contribu- 
tion from flip-flops between nuclei of the common species 
i. Assuming a uniform distribution of all isotopes in a 
2D quantum dot with a Gaussian envelope function, we 
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Appendix A: CONTINUUM LIMIT 

In this appendix, we describe how the dimensionality 
d and envelope wave function shape parameter q are de- 
fined. For further details on the definition of these quan- 
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tities, see Ref. [l5[- For a homonuclear spin system, the 
hyperfine coupling constants are given by 



A k = Av \ip(r k )\ , 



(Al) 



where A is the total hyperfine coupling constant, vq is 
the volume occupied by a single-nucleus unit cell, and 
ip(j) is the electron envelope wave function. We assume 
an isotropic electron envelope: 



ij;(r k ) =V(0)e 5 V-o, 



(A2) 



where r is the effective Bohr radius, defined as the radial 
distance enclosing N nuclear spins, and r k is the radial 
distance enclosing k spins. In d dimensions: 

d 



vol(fc spins) Vgk 



r_k 

a B 



vol ( iV spins ) Vq N 
Inserting Eqs. (|A3]) and ££2) into Eq. (|A"Ij) : 

A k = A a e~(^y /d . 



(A3) 



(A4) 



To determine the k 
normalization: 



coupling Aq, we enforce the 



Y^A k = AvoY,^)] 2 « A [ d 3 r\4>{r)\ 2 = A. (A5) 

L. 7„ J 



This gives 



A = Ac 



dke («) 



(A6) 



Making the change of variables u = (jj) q ^ d ' , we find im- 
mediately 



A = A -N / duw e 
9 Jo 



. d d 
A N-T - 

q \q 



which gives the final form for A k : 

A _i±\ q i d 

= \ N i 



A k = 



N q T (l) 



(A7) 



(A8) 



Appendix B: EFFECTIVE HAMILTONIAN 

In this appendix we give details leading to the deriva- 
tion of the effective Hamiltonian, described by Eqs. (2), 
(3), and (4) of the main text. Similar effective Hamil- 
tonians have been derived previously in Refs. [32j and 
[29l ] , but due to some differences in method and approxi- 
mation, we give additional details here for the interested 
reader. We begin from the hyperfine Hamiltonian: 



Hhf 
Ho 

V s 
h 



Ho + V s , 

{b + h z )S z +bY,l k Il 



1 



(S+h- + S-h 



+ J > 



= j2 AkIk - 



(Bl) 
(B2) 

(B3) 
(B4) 



To find an effective Hamiltonian that eliminates the flip- 
flop term Vs at leading order, we apply a unitary trans- 
formation: 



H = e s H hf e 



(B5) 



where S = —S< to ensure unitarity. We now expand 
Eq. (|B5[) in powers of S, retaining terms up to O (Vq) , 
assuming S ~ O (Vs): 

H = H + V S - [H ,S] - [V S ,S\ 

+ ±[S,[S,H ]] + O(Vi). (B6) 

To eliminate Vg at leading order, we must choose S to 
satisfy Vg — [Ho, S] = 0. The S that satisfies this relation 
is given by 



S= — V B ; L 0=[Ho,0], 

Lo 



(B7) 



which is of order Vs, justifying our previous assumption: 
S ~ 0(Vs). Re-inserting Eq. ([B7} into Eq. CE]), we 
find, up to corrections that are third-, or higher-order in 



H = H + O (V£) , 
H = H + hs,Vs] 



(B8) 
(B9) 



Directly evaluating Eq. (|B7p with Ho defined in Eq. (|B2 
and Vs defined in Eq. (|B3|) gives 



k \ 



5+j- 

b + h z + 4* - b~f k k 



1 s~I + 

b + h z -4±-b lk k 



Inserting Eq. (IBlOjl into Eq. (|B9]) gives 
H = IT) (Tl #T + U) (J-l H U 

k 

Hi = -Ub + h^ + b^lkl^-h 



(BIO) 



(BH) 
(B12) 

(B13) 



Here, the contributions resulting from the term second- 
order in Vs are given explicitly by 



h r = lj2 AkAl 



! /-/+ 

b + h* + A k /2-lrfk k 1 



b + h z - A k /2-b lk 



k > 



(B14) 



6 



or 



H = \Y. AkAl 



b + h* - A k /2-b lk Ikl1 



1+ /- 

1 b + h z + A k /2-b~/ k k 



(B15) 



We can rewrite H in terms of spin operators using 



|T)(T| = i + ^ and ll)(H 

from the main text: 



S z , which gives Eq. (2) 



H = {uj + X)S z + D, (B16) 

X = (l-P d )(fc T +/i ; ), (B17) 

£> = bj^lkl^ + l^-h), (B18) 
& 

w = b + h z + P d (h 1 +h t ). (B19) 



6 + ft 2 



4 



2 ^ b 

fe 



(/*(!* + 1) - {F k f) . (B26) 



Neglecting further corrections that are smaller by the fac- 
tor bj k /uj ~ 7fc ~ 1CP 3 in Eq. (|B23[) and terms of order 



< 



Nb 



in Eqs. (|B24[1 and (|B26|) . we arrive 



'fe b+h z —b'yk 

immediately at Eqs. (3) and (4) of the main text. The 
terms of order ~ A 2 /Nb may become important on a 
time scale r ~ Nb/A 2 . In our treatment, this time scale 
is long compared to the bath correlation time r c ~ N/A 
in the perturbative regime A/b < 1, and so neglecting 
these terms is justified. 



In the above expressions, we have introduced the diag- 
onal projection superoperator PdO = J2i 10 ('I ('I O \l), 
where the index I runs over all nuclear-spin product states 
10 = ®fc K fcm i)' ^ e now a PP^ tne commutation rela- 
tion 7 ; ~] = 2Ij:5ki and expand the prefactors in Eqs. 



(|B14|) and (|B15[) in terms of the smallness parameter 
b+iv*-b~i k ~ 77 4 ^ 1- At l ca d m g order in the expan- 
sion, we find hi t i ~ ^T°ii where 



8 Z-^ 



4^ (/-/+ 



^ = ly 



A k Ai 



8^b + h z 



bjk 



TO 



(B20) 
(B21) 



By commuting the nuclear spin operators, Eqs. (|B20I) 
and (|B21[) can be rewritten to give 



(0) 



2 ^ b + h z - b lk 



(B22) 



This relation allows us to approximate the various terms 
in Eqs. (|Bl7|l . (|B18| . and (|Hl9|> : 



(i - p d ) (2/ l ; o) 



JE 



in 



k i > 



(B23) 



Appendix C: BORN APPROXIMATION 

In this appendix we give further detail on the Born 
approximation. We begin from the equation of motion for 
the transverse spin in the rotating frame xt after applying 
the Markov approximation, neglecting the correction e(t) 
(following Eq. (7)): 



x t = — i I drS(r)a; t , 
'o 

_ — i{ui n +Aui)t\ 



(CI) 



E(t) = e- l ^" + ^ u > t 'E(t), (C2) 
E(i) = -iTrS+LQe- iLQt LQ|n) (n|5_. (C3) 

In general, it is not simple to find the exact form of the 
self energy (memory kernel) E(f). Fortunately, it is pos- 
sible to generate a systematic expansion in the perturba- 
tion V = XS Z oc 1/b, valid for sufficiently large Zccman 
splitting b > A-j^- 



E(t) = E^(t) + E< 4 )(i) + 



(C4) 



where E'"'(f) indicates a term of order ~ 0(V n ) 



O 



The expansion is performed most conve- 



niently in terms of the Laplace-transformed variable 



S(s)=£[E(i)] 



dte~ st E(t). 



(C5) 



and 



6 + 



4 



4 (6 + /i 2 -&7fe) 



2/i 



(0) 



-4? 



fe 



&7fc 



l fe> 



(B24) 



(B25) 



We expand the propagator £ [e _lLQt ] = a+ jLQ by divid- 
ing the full Liouvillian into unperturbed and perturbed 
parts: L = Lo + Ly, where Lo and Ly are defined by their 
action on an arbitrary operator O through LqO = [Ho, O] 
and LyO = [V,0]. To obtain an expansion in terms of 
the perturbation Ly, we now iterate the Dyson identity 
in Laplace space: 



HQ 



LyQ- 



zLoQ 



-0(l 2 v ). 
(C6) 
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Inserting the iterated expression (Eq. (|C6[) ) into the 
Laplace-transformed version of Eq. (|CJ3|) . we find the 
self energy in Born approximation (to second order in V) 
is 



S (2) (.s) = -iTi 



S+il- zL Q 



1 



s + iL 



xL v ———L v \n) (n\ S- 
S + iL 



(C7) 



We have simplified the above expression using the fol- 
lowing identities for the projection superoperators Q = 
1 - |n) (n\ Tr/ and P = 1 - Q: 



PL P = UP, 
PLy \n){n\ = 0, 
QL Q = QL , 



(C8) 
(C9) 
(CIO) 



which can be proven directly. To further reduce the above 
expression, we evaluate the action of L and Ly on the 
electron spin operator S- : 



LySL - --UtS"_, 



loS- 



-L+ + L 



D 



(Cll) 
(C12) 



where 5D kn = D k - D n , 5u nk — uj n -uj kl and uj k , D k are 
the eigenvalues associated with eigenstate \k): u>\k) = 
u>k\k), D\k) — D k \k). Additionally, we have denoted 
X kn = (k\X\n). 

From Eqs. flUT), ([D2]). and ([C5|, the electron-spin 
decoherence rate within a Born-Markov approximation 
will now be given by 



To 



-Im£ (2) (s = i(u n + Aw) + + ), 



(C20) 



where 



where + denotes a positive infinitesimal. Our goal here 
is to find the leading-order dependence of I/T2 on 1/6 for 
large Zeeman splitting: b > A. We therefore set Aui — 

-Re£( 2 )( S = %{u„ + Alu) + 0+) ~ O (£) 2 ) « 0, 

since this term will lead to higher-order corrections in 
1/6 within the perturbative regime. Additionally, not- 
ing that the matrix element X kn induces a flip-flop for 
spins at two sites &i,2, we find |&Dfcn| = \b{^ kl — 7fc 2 )| 
and |i5wfe n | = \A kl — A k2 \. In the case of a homonu- 
clear system 7^ = j k2 , we can set 5D kn — in Eq. 
(|C19|) . Otherwise, in a sufficiently large magnetic field 
\b (7fei - 7fc a )l > l^fei _ ^fc 2 |) we nn d a negligible contri- 
bution to the decoherence rate for terms from two differ- 
ent isotopic species (where j kl ^ 7fc 2 ), i.e., heteronuclear 
flip-flops no longer conserve energy, although homonu- 
clear fip-flops (for which j kl — j k2 ) will still occur. Re- 
stricting the sum to homonuclear flip-flops and setting 
SD nk = in this regime gives 



L+O = [X,0] + 
[D,0], 



-x 

L D 



(C13) 
(C14) 
(C15) 



and here we denote anticommutators with a '+' sub- 
script: [A, B] + = AB + BA. This leads to 



= EE I* 



(C21) 



where X 3 kr = 



(k\ X 3 



and X 3 is restricted to run over 



flip-flops between nuclei of the common species j at sites 
denoted by the indices kj , lj : 



E' 2 »(») = 



xL+- 



S + l ( L£) 



1 T7T\ L x I") ("I 



Now, noting that 



Q\n)(n\ = 0, 

Q|fc)(fc| = \k)(k\-\n)(n\, 



we can evaluate Eq. (|C16|) directly, giving 



E' 2 )( S + !Wn ) = --^|X fc „| 
/ 1 



S + iSiUnk 

1 



(C16) 



(C17) 
(C18) 



s + i (5D kn + \5uj nk ) s - i (5D kn - \8uj nk ) 



(C19) 



1 ^ A{ M 

•v 77 E ^b^Vu- 



(C22) 



Inserting Eq. (|C21[) for a homonuclear system (one iso- 
topic species j) into Eq. (|C20|) and inverting the Laplace 
transform leads directly to Eq. (10) of the main text. 



Appendix D: DECOHERENCE RATE 



Applying Eq. (jC20|) (setting Auj « 0) with Eq. (jC2T|) 
gives the rate 



j k 



which can be found directly from the formula 

}. n+ =~P-TinS(x), 
x ± iu + x 



(Dl) 



(D2) 
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where V indicates that the principle value should be 
taken in any integral over x. Rewriting Eq. (|D1|) us ~ 
ing the definition of X 3 given in Eq. (|C22j) : 



1 



IT 

~ 4 

*££ 

3 kj^l. 



c - c + 
UJkUJn 



K) a K) 9 *K 



(D3) 



where kj and are restricted to run over sites occupied 
by isotopic species j. The coefficients c± J give the ex- 
pectation value of the operator I ^ iTf; with respect to the 
initial state: 



\n) 



(D4) 



= P{P + l)-{n\F k .{P k .±l)\n). (D5) 



With small corrections of order A/Nb < 1, we can re- 
place LOk — o-> n in the denominator of Eq. (|D3|) . If the 
various nuclear isotopes are uniformly distributed with 
isotopic concentrations i/j, we allow the sum over kj, lj 
to extend over all sites fc, I at the expense of a weight 
factor Vj for each index: 

Additionally, we assume that the system is uniformly po- 
larized on the scale of variation of the hyperfine coupling 



constants so that the coefficients c^_ fe can be replaced by 
average values c?± = ((c'±\\ (double angle brackets in- 
dicate an average over all sites) and taken out of the sum. 
Finally, we change the sums over sites to a double integral 
using the prescription and coupling constants described 
in Appendix [5] neglecting the small O (1/-/V) correction 
due to the requirement k =/= I: 



dk 



dl. 



(D7) 



These approximations give 

1 _ 7T 



E"?- ; ' ■ J™ dk dl (4) 2 (4) 2 5 (4 - A\ 

(D8) 

Inserting the coupling constants defined by Eq. (|A8[) and 
evaluating the integrals gives 



T 2 4 J \qj^ 3 - 



N \uj n 



(D9) 



with the geometrical factor f(d/q) given by Eq. (13) of 
the main text. Eq. (|D9|) reduces to Eqs. (12), (14), and 
(15) of the main text in the special cases discussed there. 
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The integral in Eq. becomes undefined if the memory 
kernel has an asymptotic time dependence T,(t) ~ l/t a , 
where a < 1, and consequently the Markov approximation 
breaks down in this case. A weaker version of Markovian 
violation can occur more generally for a < 2, in which case 
the bound, Eq. @, may still be small for times t ~ T2, 
but grows unbounded in time. This situation occurs, for 
example, in the ohmic spin-boson model.— 



